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Abstract 

The observed temperature dependence of electrical conduction in polymer electrolytes is usu- 
ally fitted with two separated equations: an Arrhenius equation at low temperatures and Vogel- 
Tamman-Fulcher (VTF) at high temperatures. We report here a derivation of a single equation 
to explain the variation of electrical conduction in polymer electrolytes at all temperature ranges. 
Our single equation is in agreement with the experimental data. 



There is a wide interest in composites of ceramic nanoparticles and solid polymer electrolytes 
(SPE) for their potential applications in batteries, fuel cells, solar cells, and other electrochemical 
devices [1-4]. Usually, the SPE only exhibits a high electrical conduction at high temperatues, 
i.e., above the melting point. However, dispersing ceramic nanoparticles in the SPE matrix 
can enhance the conductivity even at low temperatures by three order of magnitudes. The role 
of the dispersed nanoparticles is to influence the recrystalization kinetics of the SPE polymer 
chain, thereby ultimately promoting localized amorphous regions and thus enhancement of the 
ion transport [2] . A number of researchs have been performed to investigate the enhancement of 
conductivity in nanoparticles/SPE composites at room temperature to a practically useful values 
(of about 1(T 4 S/cm) [2-7]. 

Understanding the temperature dependent conductvity of SPE is strictly important in order 
to realize the practically useful nanoparticles/SPE composites. Indeed, it was demonstrated that 
the general behavior of conduction in nanoparticles/SPE composites is similar to that in pure 
SPE [2-4,7]. However, a detailed understanding of the conductive mechanism is still lacking [8,9]. 

The conductivity of SPE is attributed by cations transfer. At low temperatures, where 
kT is far below the activation energy, the concept of ion conduction being the result of infre- 
quent individual ion hops over large energy barriers is an appropriate picture of the conduc- 
tion dynamics. This picture exhibits an Arrhenius temperature dependent conductivity, i.e., 
<7i(T) = (Ai/T) exp[— E a /kT\. At higher temperatures, however, the effect of temperature 
would no longer be thermally creating charge carriers for conduction, but rather to increase the 
mobility of charge carriers. This hypotesis has been used to explain the nonlinearity of con- 
ductivity in glass materials [9]. The ionic motion is strongly coupled to the polymer segmental 
relaxation to enhance the mobility. Still there is no single theory that is able to explain the 
temperature dependent electrical conduction in SPE for all temperatures. In the present work 



we develop a model that describes the electrical conduction properties in polymer electrolytes at 
all temperature ranges. 

In Fig. 1, the SPE at arbitrary temperature is divided into a number of similar cells. Each 
cell contains a number of clusters (chain segment) and each cluster is governed by a number of 
sites. Consider a certain cluster containing m sites: a part is " amorphous" sites and the rest are 
"crystalline" sites. Other authors classify as "solid" atoms and "liquid" atoms [10]. The average 
interaction energy per site can be written as E m = —(J /zm) Vi&j, with <7j = ±1 denotes 
two states: "amorphous" and "crystalline", z is the number of nearest neighbors sites, and J Q as 
a constant (interaction strength) [11]. 

Suppose there are t "amorphous" and m — £ "crystalline" sites in a cluster of m sites. 
The interaction energy per site (Bragg- Williams formalism) is e m (£) = e(L) = —J Q L 2 /2, with 
L = 2£/m — 1 and the partition function is Z m = XX { m V( m ~~ £)\£\} exp[—me m (£)/kT]. 
Using Stirling approximation and converting the summation into integral, one finds Z m = 
dL exp[-mf(L)/kT], with f(L) = e(L) + (kT/2) {(1+L) ln[(l + L)/2] + (l-L) ln[(l-L)/2]}. 
If m is sufficiently large, only the minimum free energy contributes significantly to the integration. 
This appears at L — that satisfies df /dL = or = tanh [L^Jo/kT]. For T > J a /k, the 
only solution is = 0. This standard result is only valid for the macroscopic limit [11]. For a 
mesoscopic system, the energy can be lowered as e m ~ dL(—J L 2 /2) exp[—mf(L)/kT]/Z m . 
For large but finite m the integral may be approximated by a steepest descent procedure and 
found that e m ~ — J c /2(1 — Jo/kT). The energies of clusters are similar regarless of their size. 
The energy of a cell containing n clusters can be approximated as E n = —n.J /2{\ — Jo/kT). 

The instant velocity of ion in a SPE is Vi = Vi P + v p with v ip the velocity of ion relative to 
polymer chain and v p the velocity of polymer chain. The average kinetic energy of ion become 
(KE) = m i (vf p )/2 + mi(v ip .v p ) + m i {v 2 v )/2. Since the ions are strongly coupled to polymer 



chain, the ions follow the fluctuation of polymer chain. And since the segmental motion of chain 
dominates at high temperatures, one has \v ip \ <C \v p \, so that (KE) oc (Vp), i.e., proportional to 
the kinetic energy of polymer in a cell, E n . The total energy of ion in a cell containing n cluster 
can be approximated by E l n = Vi + nE n , with k ia a real constant. The potential energy of ion, 
Vi, is dominated by interaction of cation and anion and the interaction with the nearest atom in 
polymer chain at which the cation is attached. The potential energy therefore can be considered 
to be constant, independent of temperature and chain fluctuation. 

To describe the effect of temperature on conductivity, let us start by investigating the con- 
ductance between cells. A mixture of amorphous and crystalline cells is similar to the mixture 
of conducting filler and insulator matrix. The conducting filler is associated with the amorphous 
cell and the matrix is associated with the crystalline cells. The increase in temperature resulting 
the increase in the fraction of amorphous state is similar to the increase in the volume fraction 
of conduction filler loaded in insulator matrix. 

Firstly we calculate the fraction of amorphous state. Due to fluctuation, there is a probability 
for the melting process to occur at an arbitrary temperature. Assuming a Gaussian type of 
fluctuation, the probability for melting in the range of temperature between T and T + dT is 



w(T)dT = [l/^/27r((AT) 2 ) m ] exp [-(AT) 2 /2((AT) 2 ) m ] dAT, where T m is the commonly observed 
melting point, AT = T — T m , dAT = d(T — T m ) = dT and ((AT) 2 ) m is the assembly average 
of (AT) 2 at T m [12]. The probability for the presence of the amorphous state at temperature T 
is equal to the probability for attaining melting points at all temperatures below T, i.e. p(T) = 



with erf(x) = (2/y/n) / S exp[-t 2 /2}dt. Since ((AT) 2 ) = RT 2 /C V with C v is the molar heat 



JZ» [l/v/27r((AT) 2 ) m ] exp [—(AT) 



2 /2((AT') 2 ) m ] dAT' or 




(1) 



capacity at constant temperature [12], one can approximate ((AT) 2 ) m = RT^/C vmif with C v 
is the average molar heat capacity at T m and 

f 1/2 - 1/2 erf [-C vm * AT/R] if T < T m , 



P(T) 



(2) 

[ 1/2 + 1/2 erf [C vm * AT/R] if T > T m . 



This probability denotes the fraction of amorphous phase, i.e., w a (T) = p(T) [13]. 

Based on Fig. 1, define / as the coordination number (cell hands) and a the probability for 
direct bonding of two cells. In a polycell containing n cells (n-cell), some hands formed bonds 
and the others remain free. It requires the attachment of (n — 1) new cells following a preselected 
cell to form n-cell. If only one bond is created at each attachment of a new cell, (n—1) bonds 
and (/ — 2)n + 1 free hands (excluded a freely preselected hand) appeared having the n-cell been 
formed. Since the probability of hand to form bond is a and to become free is (1 — a), the total 
probability for creating n-cell of any configuration is 

P n = n n a n -\l-a)U- 2 > +1 (3) 

with Q n is the total number of configurations. 

To form (n—1) bonds in an n-cell, it allows to select (n—1) from the total (n—1)/ hands belong 
to the added cells in [fin — l)]!/[(/ — 2)n + l]!(n — 1)! ways, and the added cells can be attached 
sequentially in (n — 1)! ways. Noting that all cells in the n-cell are identical, so that we have 
to add a division factor n! and obtaining the expression Vt n = [f(n — l)]!/[(/ — 2)n + l]!(n — 1)! 
x(n - l)!/n! = [f(n - 1)]!/ [(/ - 2)n + l]!n!. 

The alternative expression for P n is P n = Number of free hands on n-cell/The total number 
of free hands [14]. The number of free hands belongs to an n-cell (the preselected cell plus the 
added cells) is (/ — 2)n + 2. Suppose N n denotes the population of n-cell and N Q the number 
of all cells. The total number of hands belongs to all cells is N Q f and the total number of free 



hands is N Q f(l — a). Therefore 

P n = [(/ - 2)n + 2] N n /NJ(l - a) (4) 

From Eqs (3) and (4) one obtains N n = NJ[(1 - a) 2 /a]tt n f3 n , with (3 = a(l - a) f - 2 . 

The weight fraction of n-cell is w n = nN n /N D = f[(l — a) 2 / a]nVt n f3 n , with Y^=\ w n = 
Y^=i /[(l — a) 2 / a}nVL n (3 n = 1. The summation result on the right hand side depends on (3. For 
a specific value of (3, there is only one summation result. Since (3 is a polynomial function of a, 
a specific value of (3 results more than one value of a. By replacing a with a' that also satisfies 
(3 — a'(l — a')f~ 2 , the summation result on the right hand side never changes. However, since 
in general (1 — a) 2 /a ^ (1 — a') 2 /a', we have, in general, X^Li w n 7^ 1, that sounds unphysical. 
The physical soundness is achieved only when the smallest root of f3 — a(l — a)f~ 2 is used [14]. 
Therefore, the accepted expression for w n would be w n = /[(l — a') 2 / a']nVt n f3 n with a' is the 
smallest root of a'(l — a'Y~ 2 = a(l — a)^ 2 . 

The summation of fl n for all finite n gives the weight fraction of the so called sol with respect 
to the amorphous phase, w s = J2 a \i finite n w n — (1 — o) 2 ol j (1 — a') 2 a}. The weight fraction of the 
so called gel (infinity network) with respect to the amorphous phase is then 

w g = l-(l- a) 2 a'/(l - a'fa. (5) 

The weight fraction of gel with respect to the total weight of SPE is 

w' g = w a w g . (6) 

Since the presence of segmental fluctuation, the number of cluster is random in all cells 
to result the randomness in ionic energy of cells. The amorphous cells also occupy random 
positions in a SPE. Both energy and position of amorphous cells satisfy the requirement for 
the application of critical path method approach [15] so that the conductance between two cells 



having n and n' clusters can be written as g n , n i = (g /T) exp[— 2'js nn ' — W^ n ,/kT], with g Q is 
the prefactor, 7 is the screening parameter, s nn > is the distance between two cells, and W^ n , = 
k (\ne n \ + \n'e n i \ + \ne n — n'e n >\)/2. 

We only need to consider the connected amorphous cells since only them contribute to the 
electrical conduction. Since the amorphous cells are connected at a constant distance (for nearest 
neighbor ones), the first term in exponential is constant. Sheng, Sinchel, and Gittleman argued 
that the temperature dependence of conductivity of composite consists of carbon particles dis- 
persed in polyvinylchloride is similar to that of between two carbon particles in that composites 
[16]. Following this argument, the conductivity of PE (at high temperatures) is given by 

a 2 (T) = ^eM-B/k(T-T )] (7) 

with B is taken to be the average of (\n — n'\ + n + n') multiplied by kJ , and T Q = J a /k, which 
proves that at high temperatures, the conductivity varies according to VTF law [8,17]. 

Although the conductivity of sol is higher than that of crystalline phase, since it is always 
covered by crystalline phase, the effect of sol on conductivity of the PE is ignorable. It is analog 
to insulator coated conducting particles, where the presence of insulator on the surfaces blocks 
the transport of charges. Practically, the polymer can be classified into two phases: a gel phase 
with conductivity cr 2 (T) and a mixture of crystalline phase and sol with a conductivity of (Ji(T). 
The total conductivity of polymer is then 

a(T) = (l-v g )a 1 (T)+v g a 2 (T) (8) 

As an illustration let us use the above result to explain the experimental observation of 
the electrical conductivities of polyethylene oxide (PEO) based polymer electrolytes. We select 
cubical shapes of cell such that it contact with six other cell (/ = 6). The probability for direct 
bonding of two cells a is equal to w a , the volume fraction of amorphous phase. The heat capacity 



of polyethylene oxide at constant pressure is around C p = 30 J/mol K [18]. For relatively narrow 
temperature range one can assume that heat capacity only changes very little with temperature 
(can be considered to be constant). Using a rough approximation (that usually used for gas at 
low pressure), i.e., C v = C p — R, we used in the calculation C vm * = 25 J/mol K. 

Figure 2 shows the comparison of the present model and experimental data for a system of: 
(a) Cd . 7 5PS 3 Lio.5(PEO) [8], (b) Cd . 7 5PS 3 Nao.5(PEO) [8], (c) Cd . 7 5PS 3 Lio.5(PPG) [17], with 
PPG denotes polypropylene glycol. The parameters used in calculation are displayed in Table 
1. The values of those parameters are close to that used by Jeevanandam and Vasudevan to fit 
the experimental data using two separated curves: an only-Arrhenius curve for low temperatures 
conductivity and only-VTF curve at high temperatures [8,17]. Indeed, these values often depend 
on the preparation condition. Eventually, similar material provides different values of these 
parameters. The present model results in a single equation for the conductivity behavior of solid 
polymer electrolytes and also successfully explains the change in electrical conductivity from 
Arrhenius like at low temperatures to VTF like at high temperatures with a continuous slope. 

We will show that the parameters used in calculations are acceptable. By nothing that 
w a = a, it can be shown that if w a < l/(f — 1), w g = 0. According to Eq. (8), this condition 
implies that only Arrhenius behavior to occur. When w a > l/(f — 1), w g > so that the 
VTF behavior occurs. Thus the temperature at which w g changes from zero to non-zero can be 
considered as the critical temperature, at which the VTF behavior of conductivity starts to occur 
and becomes dominant by further increasing the temperature. This temperature is known as the 
glass temperature, T g , that satisfies w a {T g ) — l/(f — 1). Since the glass temperature is less than 
the melting point then using Eq. (2) one has l/(f - 1) = 1/2 - 1/2 erf \^C vm jR{T m - T g )/T m ], 



or 

~ 2 



/ /;„ ! I - J^erf 



1 - 



(9) 



R 

Using / = 6 and C vmif = 25 J/mol K, we have T g /Tm = 0.68 The ratio of T g /T rn observed on a 
lot of polymers including polyethylene, polypropylene, polystyrene, PEO are located between 0.5 
and 0.8 [19]. Since our result is located in the range of observation date, the selected parameters 
C vmif as well as / are acceptable. 

To verify the acceptability of parameters T m and T a , let us use the relation between T g and 
T Q (Vogell temperature) as T g = T Q + 50 K [17]. Using T Q in Tabel 1, the glass temperature of 
the three types polymer electrolytes are 140.7 K, 140.7 K, and 192.2 K for Cd .75PS 3 Lio.5(PEO), 
Cd .7 5 PS3Nao.5(PEO), and Cd .75PS 3 Li .5(PPG), respectively. Furthermore, using T m in the last 
column of Tabel 1, we have the ratio of T g and T m for these three SPE are 0.6544, 0.6544, and 
0.6537 for Cd . 7 5PS 3 Li .5(PEO), Cd .75PS 3 Nao.5(PEO), and Cd .75PS 3 Lio.5(PPG), respectively. 
Again, these values are also located in the experimental observation, proved that the selected T a 
and T m are acceptable. 

It is obvious from Eq. (9) that the glass temperature is located near the melting point 
if C v is very large. This behavior can be understood easily from eq. (2) and the relation 
((AT) 2 ) = RT 2 /C V , where C v is inversely proportional to the standard deviation of fluctuation, 
thereby the large value of C v implies to the small value of standard deviation such that the 
melting point is localized in a small range of temperatures. Since the factor generally depends 
on the type of polymer, the ratio of T g /T m is also polymer specific. 

In conclusion, we have shown that by considering the fluctuation of melting point and the 
application of gelation theory the determine the formation of amorphous states network we can 
explain the presence of two conductivity behavior in polymer electrolytes. We showed that the 
VTF behavior at high temperatures can be derived using a simple Ising model. Our result in 



the form of a single equation for electrical conductivity at all temperature ranges is in agreement 
with several experimental evidences. Our results are easily extended to explain the conductivity 
behavior of nanoparticles/SPE composites. 
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Figure Captions 

Figure 1. 

The modeling steps used in this work: (a) At arbitrary temperature the SPE always con- 
tains amorphous and crystalline phase, (b) The SPE is divided into a number of similar cells. 
The content of each cell is either dominated by amorphous phase or crystalline phase. We can 
consider the SPE as a "solution" of amorphous cells into a "solvent" of crystalline phase. Some 
amorphous cells form network. Networks that contain finite number of amorphous cell is named 
as "sol" and that contain infinity number of amorphous cells is named as gel. (c) The formation 
of cells network can be similarized as the polymerization process. The cell denoted by (1) is the 
preselected one and the marked cell hand denotes the preselected hand. 

Figure 2. 

The comparison of the theory calculation (lines) and the experimental data (symbols) for 
systems of (a) Cd .75PS 3 Nao.5(PEO) [8], [open circle, solid line], (b) Cd .75PS 3 Lio.5(PEO) [8], 
[plus symbol, dotted line] and (c) Cdo.75PS3L.io. 5 (PPG) [17] [cross, dashed line]. 



Table 1. The parameters used to calculate the theoretical lines in Figure 2. 
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